--- title: MAP : Module to display the data on a map keywords: fastai sidebar: home_sidebar ---
{% raw %}
{% endraw %} {% raw %} 
# Put these at the top of every notebook, to get automatic reloading and inline plotting
%reload_ext autoreload
%autoreload 2
%matplotlib inline
from IPython.core.interactiveshell import InteractiveShell
InteractiveShell.ast_node_interactivity = "all"
{% endraw %} {% raw %} 
# If the graph do not display, try :
#!jupyter labextension install jupyterlab-plotly
# You may have to upgrade Node and Jupyter
{% endraw %} {% raw %}
{% endraw %}

Display a 2D scatter plot of the mesured CO2 level together with the evaluation method

{% raw %}
{% endraw %} {% raw %}

build_sounding_scatter[source]

build_sounding_scatter(df_sounding, gaussian_param, with_dash=False)

{% endraw %} {% raw %} 
from oco2peak.datasets import Datasets
from oco2peak import find_peak
import re
import plotly.io as pio
import json

config_file = './configs/config.json'
# Load config
with open(config_file) as json_data_file:
    config = json.load(json_data_file)
datasets = Datasets(config_file)
url_month = datasets.get_files_urls('result_for_oco2_1808')[0]
url_peak = datasets.get_files_urls('peak_data-si_2018082505')[0]
print(url_peak)
df_month = datasets.get_dataframe(url_month)
#df_month.head(3)
df_peak = datasets.get_dataframe(url_peak)

https://storage.gra.cloud.ovh.net/v1/AUTH_2aaacef8e88a4ca897bb93b984bd04dd/oco2//datasets/oco-2/peaks-detected-details/peak_data-si_2018082505142073.json
{% endraw %} {% raw %} 
sounding_id = re.findall(r'(\d{16})', url_peak)[-1]
print(sounding_id)
gaussian_param = datasets.get_gaussian_param(sounding_id, df_month)
df_peak['gaussian_y'] = df_peak.distance.apply(
    lambda x: find_peak.gaussian(x=x, m=gaussian_param['slope'], b=gaussian_param['intercept'], A=gaussian_param['amplitude'], sig=gaussian_param['sigma']))

sounding_scatter = build_sounding_scatter(df_peak, gaussian_param, with_dash = False)
pio.show(sounding_scatter)

2018082505142073
{% endraw %}

Display the CO2 measurement on a map

{% raw %}
{% endraw %} {% raw %}

build_sounding_map[source]

build_sounding_map(df_sounding, mapbox_token)

{% endraw %} {% raw %} 
mapbox_token = config['mapbox_token']
sounding_scatter = build_sounding_map(df_peak, mapbox_token)
pio.show(sounding_scatter)
{% endraw %}

Display all the detected peak on a map

{% raw %}
{% endraw %} {% raw %}

build_world_map[source]

build_world_map(data)

{% endraw %} {% raw %} 
folium_map = build_world_map(df_month)
folium_map
{% endraw %} {% raw %} 
from nbdev.export import *
notebook2script()

Converted 00_core.ipynb.
Converted WIP_OCO2_Capture.ipynb.
Converted datasets.ipynb.
Converted find_peak.ipynb.
Converted index.ipynb.
Converted map.ipynb.
{% endraw %}